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Abstract 



We give an alternative method to that of Hardy- Ramanuj an- Rademacher 
to derive the leading exponential term in the asymptotic approximation to 
the partition function p(n,a), defined as the number of decompositions of a 

_ I ■ positive integer n into integer summands, with each summand appearing at 

> 

^f) • most a times in a given decomposition. The derivation involves mapping to 

an equivalent physical problem concerning the quantum entropy and energy 
ON 

currents of particles flowing in a one-dimensional (ID) channel connecting 

o' 

-I— » 

c3 



thermal reservoirs, and which obey Gentile's intermediate statistics with sta- 
tistical parameter a. The method is also applied to partitions associated with 
Haldane's fractional exclusion statistics. 



A classic result in the theory of partitions is the Hardy-Ramanujan-Rademacher formula 
for the unrestricted partition function p(n, oo), wherein the latter, combinatoric quantity is 
represented as a power series whose terms involve elementary functions of nxta This series 
yields the following asymptotic approximation: 
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p(n, oo) ~ _^WW* (1) 

A series representing p(n, 1), the number of decompositions of n into distinct summands, has 
also been derived (see, e.g., Sec. 24.2.2 of Ref. |), yielding the asymptotic approximation 

1 



pfa ^ ~ 4 . 3 i/;. W 3/4 ^ i/3 ^ ( 2 ) 

And more recentlyu Hagis used the Hardy-Ramanujan-Rademacher method to derive a 
power series representation of p(n, a) for arbitrary a = 1,2,..., yielding the asymptotic 
approximation 



p(n a) V^a 14 e VVP(i+«)l ^ ( 3 ) 

P{,) (l + a)3/4 (24n) 3 / 4 ' {) 

where n ^> a. As an example, for a = 4 the number of partitions of n — 1000 to five 
significant figures is 2.4544 x 10 28 , while approximation (^) gives 2.4527 x 10 28 , accurate to 
within 0.1%. 

In the present work, we give an alternative and more direct derivation of the asymptotic 
approximation to In p(n, a) which, from Eq. (H), is: 



lnp ^ a) ^^WTa}-^- (4) 

The derivation begins by considering a ID quantum channel which supports particles obeying 
Gentile's intermediate statistics^ characterised by statistical parameter a, the maximum 
occupation number of particles in a single particle state, with a = 1 describing fermions and 
a = oo bosons. The left end of the channel is connected to a particle source and the right 



end to a particle sink. The channel is dispersionless so that particle packets with different 
mean energies have the same velocity c and hence transmission time r = L/c, where L 
is the channel length. Imposing periodic boundary conditions on the channel length, the 
single-particle energies are €j = hfj = hj/r, j — 1, 2, . . ., where h is Planck's constant. The 
total energy E n of a given Fock state is E n = J2j e j n j — nh/r, where n = J2'jLiJ n j, an d 
rij < a is the occupation number of, say, the right-propagating mode j. 

We now suppose that the source emits a finite number of particles with fixed total energy 
E n . The maximum possible entropy of this collection of right-propagating particles subject 
to the fixed energy constraint is S(n, a) = ks In p(n, a). Thus, the problem to determine the 
asymptotic approximation to lnp(n,a) is equivalent to determining the asymptotic approx- 
imation to the entropy S(n, a) of the just-described physical system. (C.f. Sec. 4 of Ref. |7|, 
where the same set-up restricted to bosons was considered in the problem to determine the 
optimum capacity for classical information transmission down a quantum channel.) 

The crucial next step is to consider a slightly different set-up, in which the particle source 
and sink are replaced by two thermal reservoirs described by grand canonical ensembles, 
with the chemical potentials of the left and right reservoirs satisfying ^l = Hr = 0, the 
temperature of the right reservoir Tr = 0, and the temperature Tl of the left reservoir chosen 
such that the thermal- averaged energy current flowing in the channel satisfies E(Tl,o) = 
E n /r. (Note that the chemical potentials are set to zero since there is no constraint on the 
thermal-averaged particle number.) With this choice, the thermal-averaged, channel entropy 
current S(Tl, a) coincides with S(n, a)/r in the thermodynamic limit E n (equivalently n) — > 
oo. 

The advantage with using the latter, grand canonical ensemble description as opposed to 
the former, microcanonical ensemble description is the greater ease with which the energy 
and entropy currents can be calculated. The starting formula for the single channel energy 
current is: 
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E(T,a)=Y,e J [na(e j )/L]c, (5) 



where we have dropped the subscript on Tl, and where n a (e) is the intermediate statistics 
thermal-averaged occupation number of the right-moving state with energy e:u 

- n - l Q + 1 (k\ 

In the limit L — ► oo (equivalently r — > oo), we can replace the sum with an integral over 
j and, changing integration variables j — > e — (h/r)j = (hc/L)j, we have [c.f. Eq. (13) of 
Ref. |]: 
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E{T,a) = -J o deen a {t). (7) 

A formula for entropy current can be derived as follows. First note that the thermal-averaged 
occupation energy e = en a (e) and the entropy s for a given mode with energy e are related 
through the first law: ds/dT = {l/T)de/dT. Integrating with respect to temperature and 
then summing over the right propagating channel modes, we obtain 

S(T,o) = -% HdeeT d(3'(3'^. (8) 

h JO Jbeta Op' 

The integrals are straightforwardly carried out by noting from (Q) that the thermal- 
averaged occupation energy e = en a (e) of level e for statistical parameter a is just the 
difference in the thermal-averaged occupation energies of levels e and e(a + 1) for bosons. 
Thus, we require only the integrals for the bosonic case: E(T, oo) = 7r 2 (fc^T) 2 /(6/i) and 
§(T,oo)=7T 2 k%T/(3h), giving 

r/ x ( 1 \ n 2 (k B T) 2 



1 + a/ 6h 
and 
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^K 1 -!^^ (10) 

Comparing powers of T appearing in Eqs. (§) and (|10|), and recalling that E(T,a) = 
E n /r and S(T,a) ~ S(n,a)/r, we learn immediately that lnp(n,a) ~ C(a)y^n, where the 
n- independent factor C(a) is given by 



C(a) = ^ (r ' Q) . (11) 



k B \JE(T,a) 
Substituting in the expressions @ and fllQf) for _E and 5, respectively, we finally obtain 



C(a) = 7rw2a/[3(l + a)], in agreement with Eq. (H). 

We will now carry out the same steps as above for particles obeying Haldane's fractional 
exclusion statistics^ to derive the asymptotic approximation to the logarithm of yet another 
type of partition function, p(n,g), which also interpolates between the unrestricted and 
distinct partition functions [Eqs. ([I]) and (|2|), respectively]. Following the usual conventions, 
the statistics parameter is denoted by g = 1/a (so that g = describes bosons and g = 1 
fermions). Partitions associated with exclusion statistics are subject to additional constraints 
as compared with partitions associated with intermediate statistics (see below). 

The energy and entropy currents for particles obeying exclusion statistics are"'" 

E(T,g) = V^pl J™ dxxn g (x) (12) 

and 

- k 2 T f°° 

S(T, g)= f— / dx {n g In n g + (1 - gn g ) ln(l - gn g ) 

n Jo 

- [l + (l-g)n g ]\n[l + (l-g)n g ]}, (13) 

where x = (3e and the thermal-averaged occupation number isli3 

1 (x) = [w(x)+ g y\ (14) 

with the function w(x) given by the implicit equation 

w(x) 9 [l + w(x)} 1 - 9 = e x . (15) 



n 9 \ 



Again, comparing powers of T appearing in Eqs. ([12]) and (0), we learn immediately 
that lnp(n,g) ~ C(g)y/n, where the n-independent factor C(g) is given in terms of E and 
S as in Eq. fllTl) . Substituting in the expressions for E and S and performing a change of 
variables from x to u>,E3 Eq. fllTD becomes after some algebra 



O//): 4$=, (16) 
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where 

r°° 1 
e(p) = / dw—— [(1 - g) ln(l + w) + g In w] (17) 

and 

/■oo 

s(g) = dw[\n(l + w)/w-\nw/(w + l)}. (18) 

Using the identity s(p) = 2e(g), Eq. (ITS) can be further simplified to 



C(g) = ps{g). (19) 

Let us now describe some of the properties and consequences of result (|T9"P. Integral 



can be rewritten in terms of dilogarithmgla and only for certain choices of lower integration 
limit do closed-form solutions exist. For example, from fll5|) we have w g= o(0) = and 
w g= i(0) = 1 and solving the respective integrals, we obtain s(0) = n 2 /3 and s(l) = 7r 2 /6. 
Substituting these values into (|T9|) , we indeed obtain the arguments of the exponentials in 
the asymptotic approximations to the unrestricted and distinct partition functions, Eqs. (|lj) 
and (0) respectively. It is tempting to speculate that closed-form solutions to the integral 
s(g) exist only for g = 0, 1/2, 1/3, 1/4, and 1 in the interval [0, 1], since it is only for these 
rational values that Eq. ( |T5D can be solved analytically for the lower integration limit w g (0). 
For g = 1/2, we have w 1/2 (0) = (-1 + \fl)j1 and s(l/2) = vr 2 /5, so that 



lnp(n, 1/2) ~ tt^/2/5 ■ Vn. (20) 

Note that C , s= i/ 2 (= tc J2/5) < C a=2 (= 2ir/3), signalling the fact that ft(n,g) < p(n,a = 
1/g) for < g < 1, a consequence of additional constraints on the allowed partitions 



associated with Haldane's statistics. These constraints are discussed in Ref. |14[ The above, 
closed-form solutions for g = 0, g = 1/2, and 1 were obtained by solving the integral s(g) 
numerically and then noting that the result when divided by 7r 2 was rational. This method 



does not work for the g = 1/3, 1/4 cases, however, owing to the complicated form of the lower 
limits 1^1/3(0) and if 1/4(0) (they are roots of third and fourth degree polynomial equations, 
respectively). A more sophisticated method is required in order to determine whether or 
not closed-form solutions exist for these latter two cases. 
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